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ABSTRACT 

Light curves of eight BL Lac objects in the BVRI bands have been analyzed. 
All of the objects tend to be bluer when brighter. However spectral slope changes 
differ quantitatively from those of a sample of QSOs analyzed in a previous paper 
(Trevese & Vagnetti 2002) and appear consistent with a different nature of the 
optical continuum. A simple model representing the variability of a synchrotron 
component can explain the spectral changes. Constraints on a possible thermal 
accretion disk component contributing to the optical luminosity are discussed. 

Subject headings: galaxies: active - BL Lacertae objects: general - BL Lacertae 
objects: individual (3C 66A, PKS 0422+004, S5 0716+71, OJ 287, ON 231, OQ 
530, S5 1803+78, BL Lac) - quasars: general 



1. INTRODUCTION 

Multi-wavelength observations of active galactic nuclei (AGN) suggest the presence of 
various emission components with different relative weights in different AGN classes. The 
large number of parameters involved in physical models of emission mechanisms can hardly 
be constrained by single epoch observations. Multi-band variability studies are then crucial in 
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identifying and constraining individual components. Moreover, active nuclei are not spatially 
resolved in general. For this reason variability time scales, delays and correlations between 
variations in different continuum components or emission lines provide the main information 
about the location of different components. The very definition of AGN classes is in part 
based on variability properties. In fact, while quasars (QSOs) and Seyfert galaxies vary, on 
average, less than 0.5 mag on time scales of months in the optical band, blazars (i.e. BL Lac 
objects and Optically Violent Variables, OWs) may change their luminosities by more than 
one mag in a few days (Ulrich, Maraschi, and Urry 1997). Quantifying spectral variability 
will then help understanding different AGN classes by comparison. Although monitoring 
campaigns, from radio to X-ray wavelengths, exist for a limited number of objects, they 
are difficult to realize with adequate time sampling and duration. For this reason we try 
to exploit multi-band variability studies limited to the optical range, which in some case 
are available for several objects and relatively long observing periods. In a previous paper 
Trcvcsc & Vagnetti (2002) analyzed the spectral slope variability of a sample of PG QSOs 
on the basis of the light curves in the B and R bands, provided by the Wise Observatory 
group, who monitored 42 quasars for 7 years with a typical sampling interval of about one 
month (Giveon et al. 1999). A spectral variability parameter /3 = Ace/ A log /,^, f„ being the 
specific flux and a = d log f,^/ d log u the spectral slope, was defined by Trevese & Vagnetti 
(2002) who derived constraints on variability mechanisms from the average a and f3 values. 
In particular they showed that the variation of spectral slope implied by changes of the 
accretion rate are not able to explain (3 values deduced from the observations, while "hot 
spots" on the accretion disc, possibly caused by instability phenomena (Kawaguchi et al. 
1998), can easily account for the observed spectral variability. In the present study we 
analyze B, V, R, I light curves of eight BL Lac objects obtained by D'Amicis et al. (2002) in 
a monitoring campaign of about 5 years with an average sampling interval of ~ 25 days. We 
show that BL Lacs clearly differ from QSOs in their a, /5 distribution. We show that a simple 
model representing the variability of a synchrotron component can account for the observed 
a and P values and we discuss constraints which can be derived from spectral variability on 
the amount of thermal radiation contributed by a disk component in the optical band. 



2. DATA ANALYSIS 

The sample consists of 8 radio selected BL Lac objects hsted in Table 1. Observations 
were made in the period 1995-1999, in the standard Johnson- Gousins B, V, R, I bands and 
are described in D'Amicis et al. (2002). Typical photometric uncertainties range from 0.01 to 
0.03 mag r.m.s.. In table 1 the following quantities are reported: Column 1 - Name; Column 
2 - Goordinate Designation; Column 3 - Number of observations; Column 4 - average B 



-3- 



magnitude; Column 5 - redshift; Column 6 - average spectral slope; Column 7 - spectral 
variability index. For the present analysis we regard as "simultaneous" the observations in 
different bands taken during the same observing night. We considered the data dereddened 
as in D'Amicis et al. (2002). We then define for each object an "instantaneous" slope by 
a straight line fit through the B, V, R, 1 data points as a function of logu, neglecting for 
the moment any curvature of the spectral energy distribution within the optical range. A 
positive correlation of the instantaneous spectral slope a = 51og/i,/91ogz/ with flux changes 
6logf^j^ has been already discussed by D'Amicis et al. (2002). In the present paper we 
quantify spectral variations in two different ways. A first method considers, for each object, 
all the positions in the a — log f,^g diagram, f^^ is the specific flux at the effective frequency 
ub of the blue band. These points, together with the relevant regression lines are shown in 
Figure 1 for each of the 8 objects. In the same figure the solid lines represent the average 
increase of a with log f,^g . Clearly this quantity is independent of the time at which a and 
log fug values are measured. We refer to the blue band, instead of the red one used by 
D'Amicis et al. (2002) in their correlation analysis, for comparison with the study of QSO 
spectral variations by Trevese & Vagnetti (2002). The slopes b of the regression lines of each 
object can be assumed as a first spectral variability index. 

A second method, adopted by Trevese & Vagnetti (2002) considers, for each object, the 
quantity: 

which represents the change of spectral slope per unit variation of log fiux as a function 
of the rest frame time lag Tjj and again we refer to the fiux changes in the blue band. 
We stress that /? is a very sensitive indicator of spectral changes but, for this reason, it 
is also a noisy quantity. For flux changes comparable or smaller than the photometric 
noise, (5 values are randomly spread over a wide range. Thus we limit the computation to 
\\og fug{tj) — log fug{ti)\ > 0.03. In any case, instead of individual values, it's worth 
referring to average values, e.g. in bins of time lag r. 

To constrain physical models of variability, it's interesting to see whether spectral 
changes are different on different time scales. Figure 2 shows for each object Pij as a function 
of Tij. The mean values of Pij in bins of 200 days are also shown. The errors represent the 
standard deviation of the mean, which has been computed according to Edelson & Krolik 
(1988), to take into account correlation among Pij values within each bin. For each object 
the value of the first spectral variabifity index b is also shown for comparison. In general the 
larger variations of P occur at larger time lags where, due to the worst statistic, they are less 
reliable. For r < 500 d there are no evident systematic trends of /? as a function of r. Thus 
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it makes sense to define a less noisy index taking the average of /3 within the first 500 days: 

13= {13,,), r,, < 500 d. (2) 
In the following we limit our analysis to the average variability index (3 . 

Besides the spectral slope a computed by a straight hne fit through the log /^,. versus 
log i/j (i=B,V,R,I) data points, we also computed two other slope values, towards the blue 
and red side of the spectral range: 

0.4[(^ -R)- {Bo - K)] o.4[(y -I)- (K - Ro)] . 

l0g(AB/AR) l0g(Ay/A/) 

where Bo,Vo, Ro, lo and As, Ay, Ar, A/ are the zero points and the effective wavelengths of 
the four photometric bands (Cox 2000). For all but two objects the values of ayi and ubr 
are respectively slightly larger and slightly smaller, on average, than a, indicating a little 
curvature of the spectrum, and the slopes of the linear regressions asR ys. a and ayi vs. 
a do not differ significantly from one. In the case of BL Lac and OQ 530, ayi is almost 
equal to aBR, or slightly smaller, and becomes significantly smaller than a (and when 
a is small. A possible explanation is that the (constant) contribution of the host galaxy 
emerges in the faint phases causing a more marked softening of the spectrum on the red side 
of the wavelength interval considered. This is consistent with the fact that the host galaxy 
contribution, as reported by Urry et al. (2000), is relatively strong for these two objects. 



2.1. Compcirison with QSOs 

The average values (3 of the variability index, for all time lags smaller than 500 days 
rest-frame, are plotted in Figure 3 versus the relevant average a values. For comparison the 
values for the sample of 42 quasars discussed by Trevese Sz Vagnetti (2002) are also shown^. 
BL Lacs and quasars are clearly segregated in the a — (3 plane. BL Lacs show, on average, 
lower a values, consistently with the absence, or lower relative weight, of the thermal blue 
bump component. The new information in Figure 3 is that BL Lacs have (3 values smaller, 
on average, than quasars. BL Lac variability is stronger, a changes are also larger, but 
their ratio j3 is smaller. Moreover, using both a and j3 parameters together allows a clear 
distinction of the two classes. We want to discuss whether this behavior fits into the currently 
accepted models describing the different emission mechanisms of blazars and QSOs. 



values for QSOs are slightly different from those shown in Trevese & Vagnetti (2002) since they have 
been recomputed at a fixed rest-frame time lag interval, and taking the average over 500 instead of 1000 
days for consistency. 
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3. VARIABILITY OF THE SYNCHROTRON COMPONENT 

3.1. A simple model 

The continuum spectral energy distribution of blazars from radio frequencies to X and 
7-rays can be explained by a synchrotron emission plus inverse Compton scattering (Ghis- 
ellini & Maraschi 1992; Sikora, Begelman, & Rees 1994). Variability can be produced by an 
intermittent channeling into the jet of the energy produced by the central engine. Spada 
et al. (2001) have considered a detailed model where crossing of different shells of material, 
ejected with different velocities, produce shocks which heat the electrons responsible for the 
synchrotron emission. The resulting spectra are compared with multi-band, multi-epoch 
observations of 3C 279 from radio to 7 frequencies, showing a good agreement. Similar com- 
putations have been done for Mkn 421 (Guetta et al. 2002). In the case of the eight objects 
of our sample, B, V, R, I bands are sampling variability of the synchrotron component. This 
component can be roughly described by a broken power law characterized by the break (or 
peak in vL,^) frequency i/p and the asymptotic spectral slopes ai and 0:2 at low and high 
frequency respectively, cxi > —1, a2 < —1 (Tavecchio, Maraschi, & Ghisellini 1998). An 
equivalent representation is: 

L^^Lp- 2[(-)-"i + (-)-"2]-i = Lp ■ H{u/up- ai, 0C2) (3) 

where Lp is the specific luminosity at Vp. We adopt this representation for a stationary com- 
ponent of the SED and we add to it a second (variable) component with the same analytical 
form but different peak frequency Up/ and amplitude Lp> to produce spectral changes: 

Lp - H{v/vp; CKi, + Lpi ■ H{v/vpi] cti, (4) 

The addition of the second component mimics the behavior of the synchrotron emission 
of model spectra when shell crossing occurs, producing an increment of emission with Vpi > 
Up, due to newly accelerated electrons. With such a representation we can compute a and 
/3 as a function of Vp'/i'p, for different values of Pp and given values of ai, a2 and Lpi/Lp. 
We adopted typical values of the asymptotic slopes ai — —0.5, a2 — —2, we assigned to Up 
different values in the range 2.5-10 -10^^ Hz, and we made Vp'/vp vary in the range 0.8-4. 
We also adopted a ratio Lpi/ Lp corresponding to a 0.5 mag change in the blue band, which 
is representative of the r.m.s. variability. The results are shown in Figure 4, where we 
use as a the slope of the stationary component. The three different lines are computed for 
logz/p = 14.4,14.7,15.0 from left to right respectively. The computed curves fall naturally 
in the region occupied by the data, namely it is possible to account for the position of the 
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objects in the a — (3 plane with typical values of ai, q;2, ^'p'/i'p and Up, required to represent 
the overall SED of the objects considered (see Fossati et al. 1998). 

3.2. The Thermal Bump 

BL Lac objects can be classified according to the break frequency of synchrotron emis- 
sion, from Low frequency peaked BL Lacs (LBL) which show a break in the IR/optical 
band to High frequency peaked BL Lacs (HBL) with the synchrotron break at UV/X-ray (or 
higher) energies (e.g. Padovani & Giommi 1995). On average the blue bump (TB) produced 
by the accretion disk tends to be negligible in HBLs while in LBLs its contribution may 
be important. According to Cavaliere & D'Elia (2001) this can be interpreted in terms of 
higher accretion rate rh in the case of LBLs, producing a stronger thermal emission from 
the disk, while HBLs are dominated by non-thermal energy, partly extracted from rotational 
energy of the central Kerr hole. 

In particular the TB remains negligible during changes of the synchrotron component in 
the case of Mkn 421, a classic HBL, while emerges clearly during the faint phases in the case 
of 3C 279, a low frequency peaked object (see Spada et al. 2001; Guetta et al. 2002). We 
stress that the TB is relatively "narrow" , since its emission is mostly concentrated in the UV 
region, while the synchrotron emission extends from radio frequencies to the UV (or X-rays). 
Then the change of relative intensity of the two components can produce strong changes of 
the spectral slope in the optical- UV region, where the synchrotron and TB components have 
very different slopes (~ —2 and ~ -1-1.5 respectively). As a result, in the case of LBLs 
we could find strongly negative j3 values in the faint phases, namely a hardening instead 
of a softening of the spectrum, at wavelengths smaller but close to the peak of the TB, 
as observed in the case of 3C 279, where we can estimate /3 ~ — 1 at A ~ 2500 A on 
the basis of the data reported by Plan et al. (1999) (see their Figure 4). However these 
negative (5 values do not seem to be typical in the case of our sample, although synchrotron 
breaks are distributed in the frequency range 10^^ — 10^^ Hz. Thus, in principle, we can 
set upper hmits to contribution of the accretion disk to the overall SED on the basis of the 
observed (3 values. For a rough estimate it is sufficient to approximate the TB with a black 
body emission. As already noted in §1, although the TB itself may vary, possibly due to 
instabilities which produce hot spots in the accretion disk (Trevese & Vagnetti 2002), the 
amplitudes of its variations arc smaller than those of blazars (e.g Trevese et al. 1994). Thus 
we can approximately consider TB emission as a constant during the much larger variations 
of the synchrotron component. 

Thus we consider a specific luminosity Lg = L{u,t) of the synchrotron component with 



- 7- 



a time dependent local slope as{t) at the observing frequency 1/ — P, and a constant thermal 
TB luminosity L^b with local slope axs- The total specific luminosity at two times t and 
to will he L = Lg + Ltb and Lg — Lso + Ltb respectively. The local spectral slope a of the 
composite SED is: 

dlogL Ls Ltb ,r-^ 

a = — = —as H —aTB (5) 

a log!/ L L 

Defining y = e = and 5 — as — aso, where the subscript "o" indicates the initial 
time, we have fis — 5 /logs for the "unperturbed" synchrotron component and: 



(3= {a- ao)/\og^ 

J-in 



+ 7 — \ — VTT"^ — \y^o ~ °^TB) 



^ + x (£ + x)(i + x) 



for the composite spectrum. In Figure 5 wc report the spectral variability index /3 as a 
function of x for different values of the temperature T. 

We assume a typical value ago — —1.75 and we assume for the disk a black body 
temperature T — lO^X or T = 2 x lO^X. In both cases we compute /3 as a function of x, 
for two values of e and a typical value /S^ = 1. The result is shown in Figure 5. For 3C 
279 we can estimate (from Pian et al. 1999) ~ — 1, % ~ 1 and T ^ 2 ^ lO'^K which is 
consistent with the relevant curve in the figure. It appears that e.g. the condition /3 > 0.5 
implies x ^ 0-4 for T = lO^i^, or x < 0.2 for T = 2 x lO'^A'. Similarly (3 > implies 
X < 1 for T = lO^i^, or X < 0.5 for T = 2 x IQ'^K. For a given z/p, the ratio between 
TB and synchrotron components may, in fact, vary from typical BL Lacs to flat spectrum 
radio quasars (FSRQ), which also show intense broad emission fines. This different behavior 
depends both on the intrinsic luminosity of these components and on the viewing angle. 
An estimate of the isotropic emission can be deduced from the broad line region luminosity 
Lblr (Celotti et al. 1997) and compared with the kinetic jet power (D'Efia et al. 2003), 
as done by Padovani et al. (2003). This information is crucial in studying the relation 
among different classes of objects and testing possible AGN "grand unified theories". In this 
context, spectral variability indexes, possibly measured for stastistical samples at different 
time lags, luminosities, and rest-frame frequencies, will provide independent information on 
the relative importance of isotropic and beamed components. 



4. SUMMARY AND CONCLUSIONS 



We have performed a new analysis of the light curves of eight BL Lac objects monitored 
in B,V,R,I bands for a total period of about 5 years with an average sampling interval of 
25 days. 
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Adopting a spectral variability index /3, representing the average ratio between the 
change of spectral slope and the logarithmic luminosity change, we find that BL Lacs show 
smaller spectral variability with respect to quasars (analyzed in a previous paper). This 
happens despite the fact that both the spectral slope and luminosity variations of BL Lacs 
are larger then those of QSOs. BL Lacs and QSOs appear clearly segregated in the a — f3 
plane. Our analysis allows a quantitative comparison of the observations with a model based 
on variability of the synchrotron component. The model easily accounts for the observed 
position in the a — (5 plane with a natural choice of the relevant parameters. Thus the 
segregation in the a — (5 plane is a consequence of the different emission mechanism in the 
optical band: synchrotron in the case of BL Lacs and thermal hot spots on the accretion disk 
in the case of QSOs (see Trevese & Vagnetti 2002). A strong thermal bump can produce 
negative (3 values, which are not observed in our sample, while they do in fact appear in some 
typical LBL object, where the thermal bump emerges in the faint phases. Thus spectral 
variability, even restricted to the optical band, can be used to set limits on the relative 
contribution of the synchrotron component and the thermal component to the overall SED. 
Further studies of the spectral variability index /3 in different objects will allow us to estabhsh 
whether the thermal bump is significant only in the most extreme LBL, or to what extent 
its contribution, possibly correlated with the accretion rate m, is related to total luminosity. 

We wish to thank Alfonso Cavaliere and Enrico Massaro for discussions and comments, 
and Silvia Sclavi for providing details on photometric calibrations. This work was partly 
supported by MIUR under grant COFIN 2001/028773. 
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Fig. 1. — Instantaneous spectral slope ot versus specific flux at the effective frequency of the 
observed blue band, for the eight objects of the sample. The regression lines are also shown. 
The slopes 6, reported in the relevant panels, can be adopted as a first spectral variability 
indicator. 
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Fig. 2. — as a function of Tij for each object in the sample. Large dots represent the 
mean values of f3ij in bins of 200 days. Dashed lines correspond the first spectral variability 
indicator b. The values of b and (3 are also indicated. 
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Fig. 3. — The average values (3 of the variabihty index versus the average spectral slope 
for the eight objects of the BL Lacs sample {open circles) and the 42 QSOs (see Trevese & 
Vagnetti 2002) {dots). 
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Fig. 4. — Spectral variability index (3 versus spectral slope. Open circles represent BL Lac 
objects as in Figure 3. Lines represent the models. Each line corresponds to a fixed value 
of Vp-. \ogVp = 14.4,14.7,15.0. Along each line the ratio Vpi/vp varies from 0.8 to 4 (from 
bottom to top). We adopted an Lpi / Lp ratio corresponding a magnitude change of 0.5 mag. 
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Fig. 5. — Spectral variability versus the ratio between the thermal bump and synchrotron 
components at the observing frequency: y = The curves have been computed for two 

J^so 

values of the thermal component temperature T = 2 ■ 10^ K {thick lines) and T = 10^ K 
[thin lines), and for e = 1.1 {continuous lines) e = 1.5 {dashed lines). 
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Table 1. The sample of BL Lac objects. 



Name 


Coord. Dcsig. 


Nobs 


B 


z 


a 


/3 


3C 66A 


0219+428 


66 


15.01 


0.444 


-1.08 


0.63 


PKS 0422+004 


0422+004 


44 


15.00 




-1.05 


0.48 


85 0716+71 


0716+714 


118 


14.51 




-1.14 


0.21 


OJ 287 


0851+202 


44 


15.83 


0.306 


-1.52 


0.19 


ON 231 


1219+285 


101 


14.50 


0.102 


-1.10 


0.68 


OQ 530 


1418+546 


74 


16.22 


0.152 


-1.50 


1.03 


S5 1803+78 


1803+784 


84 


16.07 


0.680 


-1.61 


0.17 


BL Lac 


2200+420 


120 


15.56 


0.069 


-1.70 


-0.04 



